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Аннотация. Указан новьій способ линеаризации нелинейной задачи гравимет- 
рий для определения контакта. Определено, что известнье аналитические конструкции 
для определения плотностньх контактов являются его частньми случаями. Указан ряд 
поправок для последовательного итерационного уточнения аналитических аппроксима- 
ций контакта. Даньт оценки мерьт уточнения контактов в сравнений с известньюми. 
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Вступ. В процессе становления новой парадигмью в теорий интерпретации 
потенциальньх полей |5| актуальна разработка нового математического аппарата для 
моделирования потенциальньх геофизических полей, адекватного современной геофи- 
зической практике. Об адекватньїжх требованиях к новому математическому диалекту 
геофизического моделирования и ключевой роли аналитических аппроксимаций геоло- 
гической средьт и геофизических полей много сказано в работах Страхова, например в 
6). Пополняя имеющиеся наработки в зтой области, предлагаєм новьге аналитическиє 
аппроксимации для горизонтально-слоистой геологической средь. Они исследовань! в 
рамках подкласса известного класса контактньх поверхностей Страхова (31. 


О постановке задачи. Решениє 0 (х) нелинейной обратной задачи потенциала 
для контактной границь на классе Нумерова Миа) (, п) границ 2 (х) раздела одно- 


родньх сред исследовано еще в |7|. С вьсокой точностью его можно заменить на 
линейное уравнение, которое описьваєт аналитическоє продолжение вертикальной 


составляющей напряженности поля и(х) силь тяжести в сторону тяготеющих масес на 


среднюю глубину Пас (х) ІЗ: 
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Анализируя физическую суть процедурьг (1) при ограничениях, свойственньх классу 
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Вв;ражение (2) приближаєт решение задачи с точностью до квадрата разниць значений 





сі -Ю іх 
поля Й ор ) - на уровнях у-0 й у--Лп. Но зто вьгражение не позволяєт 
определить поле их п), просто продлив и(х) на уровень у-/»0. Ведь прямая у-/П 
рассекаєт тяготеющие массьт на две части, и одна из частей при 0 (х) є п будет ниже 


уровня у - п. Вследствие зтого, функция и(х, у) при у «п - негармоническая. 
О решении задачи. Найти приближениє С (п) можно не по формуле (2), а из 
предела последовательности г й (хп)) решений, полученньїжх из итераций способа 
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с последующей регуляризацией по Лаврентьеву (что достаточно для первого прибли- 
жения). Доказано |7|, что в классе ми, ) решение нелинейной контактной задачи 


С (х) с точностью П" (ре б кУ при п -» оо получаєется из линеаризованного вьтражения 


селу пн По ржехМо а) 


Из зтой формуль, при п 0,1,2 как частньве случай, получаєм известньєе анали- 


тические формульт Нумерова |2), Хьюза |8| 1 Страхова (4), где ди(2) з и(20)-- и(х.0): 
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При зтом формульт (4) и (6), (7) при п23 незффективнь для получения при- 
ближений контакта - через нарастание с номером итерации погрешностей округлений, 
особенно при вьтчислениях в окрестности точки 2 - х. Но подобньй зффект прогрес- 


-о 


сий погрешностей не свойственньй вьтражению (5). 

Во избежаниє нежелательного накопления погрешностей идбем в обход - ядра 
интегралов суммируем с разньми знаками в направленими возрастания биномиальньх 
козффициентов. Практика вьтчислений на простьїх численньїх моделях показала, что 
суммирование злементарньтх дробей в ядрах конструкций с индексами п» 0 по схеме 
(4) даєт положительнье ядра преобразований только при п- 1,23, апри п-4 - уже 
отрицательньге, что указьтвваєт на расходимость соответствующих интегральньгх квад- 
ратур. 

Итак, по результатам численного моделирования определено, что для практики 
вьічислений приближений (4) контакта 0 оп) достаточно, взяв "нулевьім прибли- 


жением" вьтражение (5), ограничиться одной из формул: 


пон? чЕ- хУВБ 
С9(хп)ви 9-2 НЕ ХУ Каву «(2-х -ли(2ИЄ, 


отця. 1Ї 86/13 хо) (2-х) р 
собача Пе і мамо в 





(2-х рлу «(2-х ау нах 





Вьбор детальности определения приближения зависит от характера исходного 


поля и(х) и масштаба построений (мерьг задачи: локальная либо региональная). Отме- 
тим особо справедливость указанньтх вьше конструкций только для областей малой ме- 
рьг («2 100 км') с несложной геометригей. 

Уточнение модели контакта. Приближение (4) 0 он) можно существенно 
уточнить с помощью ряда поправок, каждая их которьгїх есть аналитической аппрок- 
симацией искомого контакта с различной мерой детальности на классе ми (1, П) 
допустимьх искомьгх решений. 

Например, поправка типа ; | ,п)а С ох, п) ДС (п), где 


пУАРАЛЕ о це,п)- оо Іе, (9) 


ПОВЬІСИТ ТОЧНОСТЬ приближений до п? (2). Но численное моделирование показало, 
что в условиях слабоградиентного поля достаточно использовать поправку попроще 


АС (ж п)- ц(е,) Р я (10) 


где и(2,) - некое "среднее" (осреднеінное фильтрацией) значение поля на поверхности 
наблюдений. 





Другая поправка 20, РЯД щ У Ме , (оп) жд (х) , Кз 0,1,2,3.., для уточнений 
приближения С п) контакта получаєтся из решения линейного интегрального урав- 
нения 1-го рода 
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Поправку (11) можно обобщить в виде конструкции 2, (х) ша, (х)чь Ле, (х), где 


де. (х) те йт Ас, т) получена как предел последовательньх итераций вида: 
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В итоге, получен ряд аналитических конструкций (8-12), применение которькх к 
известньтм аппроксимационньм изображениям контакта (4-7) на порядок повьшиаєт 





точность их определения. Доказано, что в классе Ми) ( п) последовательнье приб- 
лижения (3-9) и (11-12) единственньт и устойчивьт (сходятся). Численно каждьтй способ 
сводится к решению линейного интегрального уравнения Фредгольма 1-го рода с яд- 
ром типа Пуассона. Анализ практической точности зтих конструкций в сравненимй с 
известньтми методами |1| возможен после апробации на полевьгх данньхх. 


СПИСОК ЛИТЕРАТУРЬ 


1. Гравиразведка. Справочник геофизика / Под ред. Е.А. Мудрецовой. - Москва: Недра, 
1981.-400 с. 

2. Нумеров Б.В. Интерпретация гравитационньх наблюдений в случає одной контакт- 
ной поверхности // ДАН СССР. - 1930. - Мо21. - С. 569-574. 

3. Старостенко В.Й., Черная Н.Н., Черньй А.В. Условия существования решения об- 
ратной задачи теории потенциала для контактной поверхности // Докл. АН УССР. 
Сер. Б. - 1988. - Моб. - С.30-33 

4. Страхов В.Н. Об интегральньх и функциональньїх уравнениях некоторьх обрат- 
ньх задач теории логарифмического потенциала и их значений для интерпретациий 
гравитационньжх и магнитньїх аномалий // Изв. АН СССР. Сер. Физика Земли. - 
1976. - Мо3. - С. 54-66. 

5. Страхов В.Н. Смена парадигмь в теории линейньх некорректньх задач. - Москва, 
2001. 

6. Страхов В.Н. Об зффективньгх по бьітродействию и точности методах построения 
линейньх аналитических аппроксимаций в геофизике, геоинформатике и гравимет- 
рим // Геофиз. журн. - 2007. - 29, Хо 1. 

7. Черньшй А.В, Дубовенко Ю.Й. Уточнение некоторьх способов приближенного опре- 
деления контактной границьт // Докладь  НАН Украийнь. - 2002. - Хо 12.- С. 128- 
132. (укр.). 

8. Ниєйез Р. Тре апайіс Ба5із ої єгаміїу іпіегргегацоп // Сеорбубз. - 1942. - 7, Хо 2. - Р. 
169-178. 








ОМ 5ОМЕ АХАГТУТІСАЇІ, АРРКОХІМАТІОХ5 КОК ТНЕ 
АРРКОХІМАТЕ РЕКІХІТІОХК Ок ТНЕ рЕХ5ІТУ СОХТАСТ5 


Рибоуепко Уигії 
Іл5цішіе ої Сеорпузіс5 МА5О, Куїу 
етаї!: петізраусуе ФикКг пеї 
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